Evaluation of Reverse Monte Carlo Models based on Molecular Dynamics 
Simulations: A Case Study of Ion Conducting Network Glasses 
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We investigate the quaUty of structural models generated by the Reverse Monte Carlo (RMC) 
method in a typical application to amorphous systems. To this end we calculate surrogate diffraction 
data from a Li20-Si02 molecular dynamics (MD) simulation and use the total scattering function, in 
addition to minimal pair distances and coordination numbers of silicon (oxygen) to oxygen (silicon) 
ions, as input for the RMC modeling. Then we compare partial radial distribution functions, 
coordination numbers, bond angles, and ring sizes predicted by the RMC models with those of the 
MD system. It is found that partial distributions functions and properties on small lengths scales, as 
distributions of coordination numbers and bond angles, are well reproduced by the RMC modeling. 
Properties in the medium-range order regime are, however, not well captured, as is demonstrated 
by comparison of ring size distributions. Due care therefore has to be exercised when extracting 
structural features from RMC models in this medium-range order regime. In particular we show 
that the occurrence of such features can be a mere consequence of the chosen starting configuration. 

PACS numbers: 61.43.Bn,61.43.-j,61.43.Fs 
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I. INTRODUCTION 

The Reverse Monte Carlo (RMC) method is commonly 
used to build structure models based on experimental 
data. Introduced by McGreevy and Pusztai in 1988,^ it 
has been spreading fast and is now considered a stan- 
dard method in analyzing structural data. Advantages 
of this method are its easy implementation and its wide 
applicability. It has been used to model various ma- 
terial systems such as crystals, polymers and glasses. 
In principal any structural data can be used as input 
for the RMC method, but most modelings focus on us- 
ing diffraction data obtained from X-ray and/or neu- 
tron scattering. In ion conducting glass systems RMC 
models of the structure of the structure have been cre- 
ated for xU2S+{l - x) SiOa glasses,^ O.TSiOa+O.SNaaO 
glass^ xNa2S-H(l -a;)B2S3 glasses^ and 0.5Li2S+0.5[(l - 
a;)GeS2+xGe02] glasses^ among others. 

As pointed out by McGreevy, RMC models are "nei- 
ther unique nor 'correct' " , but can aid our understanding 
of local structure properties and their relation to other 
physical properties.*' Accordingly it is important to qual- 
ify RMC models for different material classes and to get 
insight into the limits of this method. This is becoming 
a more urgent question now, since in recent RMC studies 
not only the short range order of various network glasses 
has been investigated, but also the medium range order. 
Among those studies are discussions of the rings sizes 
in vitreous Si02 and Ge02,^ a detailed investigations of 
amorphous GeSe2,^ and a proposal of a structural model 
for multi-component borosilicate glasses, where partial 
segregation of silicon and boron rich regions is predicted.^ 
It was also suggested to use such models as basis for fur- 
ther investigation of possible conduction pathways of the 
mobile ions. In this respect the RMC models have been 



employed in connection with geometric constraints and 
the Bond Valence (BV) analysisJ^iii (see Ref. [l^ for a 
critical discussion of this procedure) . 

In this paper we test the RMC method against struc- 
tural data obtained from a Molecular Dynamics (MD) 
simulation of a Li2 0-Si02 glass. For this purpose we cal- 
culate surrogate diffraction data from the simulated MD 
structures and these surrogate data are used as input for 
the RMC modeling. For the evaluation of the resulting 
RMC models we determine how well various properties of 
the original MD structure are reproduced. Particularly 
we compare properties such as partial radial distribution 
functions and ring-size distributions, which are not eas- 
ily accessible by experiment. Through our evaluation it 
can be clarified how far one can use the RMC method to 
gain insight into these properties, and where one has to 
be cautious to take features of the RMC model for real. 

We want to stress that for the testing to be valid, it 
is not necessary that the MD simulation is a particu- 
larly good representation of the real lithium silicate glass. 
Rather, the MD structure can be seen as a valid glass sys- 
tem in itself. We will show that the RMC models gener- 
ally compare well with the MD structure, but that one 
has to take care when analyzing features of the medium- 
range order regime. 



II. 



MOLECULAR DYNAMICS SIMULATIONS 



We perform MD simulations of a lithium silicate glass 
with the chemical formula Li2 0-Si02 using the potential 
model of J. Habasaki and I. OkadaJ^ The cubic sim- 
ulation box has a length L — 50.04 A and contains 
11664 atoms (3888 Li, 1944 Si, 5832 O) corresponding 
to a density of 2.27 g/cm'^ and a number density of 
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TABLE I: Potential parameters for the MD-simulations (cf. 
Eq. (m). 



Ion 


z 


a [A] 


MA] 


c [AVkJ/mol] 


Li+ 


0.87 


1.0155 


0.07321 


22.24 


Si^+ 


2.40 


0.8688 


0.03285 


47.43 




-1.38 


2.0474 


0.17566 


143.98 




,/o = 


4.186 kjA" 


Vol-i 


rc = 1.3A 



po = 0.093 A~^. Periodic boundary conditions are used. 
The simulations are performed in the NVE ensemble (mi- 
cro canonical ensemble where the number N of particles, 
the volume V of the simulation box, and the total energy 
E are kept constant). The energy E was adjusted so that 
the temperature of the system fluctuates around a mean 
value of 301 K with deviations of 2 K. The systems are 
equilibrated for about 1 ns and the runs for obtaining 
data have a duration of 2 ns, using a time step interval 
of At = 1 fs. 

The effective interatomic interactions between two 
atoms of type i and j at distance r are: 



■+/o(&i + ^j)exp 
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(1) 

where the parameters listed in Table U have been 
optimized— and shown to give good agreement with ex- 
perimental dataF^ii^ii^ii^. The interaction potential in 
Eq. Il]) is composed of three terms. The first one in ([1]) 
is the Coulomb interaction with effective charge num- 
bers for the species. The second term is a Born-Meyer 
type potential, which takes the short-range repulsive in- 
teractions into account, and the third is a dispersive van- 
der-Waals interaction. It is only used for interactions 
involving oxygen. 

The system was prepared by putting the atoms on a 
cubic crystal lattice and assigning to every atom random 
velocities drawn from a Maxwell-Boltzmann distribution 
corresponding to a temperature of 2500 K, which is well 
above the (computer) glass transition temperature of this 
system. From this liquid state the system is cooled down 
in several steps with intermediate periods of equilibra- 
tion. First an NVT run (canonical ensemble, where the 
number TV of particles, the volume Y and the tempera- 
ture T are fixed) of 10 ps at 2500 K is performed, fol- 
lowed by an NVE run of the same duration. After sim- 
ulating another 20 ps in the NVT ensemble and 10 ps 
under NVE conditions the temperature is decreased in 
four subsequent sequences down to 300 K. Each cooling 
cycle consists of a 10 ps run using a thermostat to de- 
crease the temperature linearly, a 10 ps NVT run at the 
target temperature, and a 10 ps NVE run to verify that 
there are no temperature drifts. The configurations at 
the end of the 300 K cooling cycle are used as starting 
points for a 800 ps long equilibration run using the NVE 



TABLE 11: Average bound coherent scattering lengths and 
atomic form factor for Li, Si, and O. 





Li 


Si 





h 


-1.9 


4.1491 


5.803 


f 


3.005 


14.41 


8.144 



ensemble. The measuring runs are 2 ns long. All MD- 
simulations were carried out with the LAMMPS software 
package^. 

The partial and total radial distribution functions 
gij{r) and G'^^^(r) as well as the total scattering struc- 
ture factor S(Q) were calculated according to the PDF- 
formalism (see [l^ for a discussion of different possible 
definitions of scattering functions). The partial radial 
distribution functions are given as 



(2) 



where (r) is the average number of particles of type j 
between distances r — dr/2 and r -f dr/2 from a particle 
of type i, and pj is the mean number density of particles 
of type j. 

The total radial distribution function is calculated by 



G 



PDF 



(wijg.j) - 1 



(3) 



where po is the total number density of the system, m 
is the number of particle types, and Wij are weighting 
factors: 



Wi 



Cihi j ^ CiCjbibj 

\i=l / = l 



(4) 



Here Ct = Pi/po are the molar fractions of particles of 
type i, and bi is their average bound coherent scattering 
length. In order to calculate X-ray diffraction functions 
one has to replace the bi with the atomic form factors fi. 

Finally, the total structure factor S{Q) is calculated 
from GP°P(r) by 



Q[SiQ) - 1] = / G^""(r) sin Qrdr 



(5) 



All data from the MD system was averaged over 11 
configurations from the 2 ns measurement run, which are 
200 ps apart each. The scattering lengths and atomic 
form factors used in Eq. 0] were taken from Ref. 19 and 
Ref. 12^, respectively, and are listed in Table |TT1 In the 
following we take the freedom to speak about these sur- 
rogate diffraction data simply as "diffraction data" and 
ask the reader to keep in mind that the data was not 
measured but calculated from Eqs.[2]to[51 
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III. REVERSE MONTE CARLO MODELING 
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FIG. 1: (color online) Comparison of G^°^(r) of Li20-Si02 
RMC models based on crystalline (black dashed line) and 
random (blue dotted line) starting configurations with that of 
a Li20-Si02 MD system (red solid line). Neutron diffraction 
data is shown in panel (a) and X-ray diffraction data in panel 
(b). 



RMC simulations were carried out using the RMC+-I- 
packagei^- We started with building an initial atomic 
configuration and then refined it against the S{Q) and 
G^^^ {r) data computed from the MD structure. 

Two starting configurations were considered: a "ran- 
dom distribution" of atoms^i (subsequently referenced 
by "R") and another from the Li20-Si02 orthorhombic 
crystal structure (space group Cmc21, subsequently ref- 
erenced by "C"). Both configurations consists of 3000 
atoms (1000 Li, 500 Si, 1500 O). The system was cho- 
sen to be cubic with a side length of 31.82 A so that the 
atomic number density is the same as in the MD struc- 
ture. From the random and crystalline starting config- 
urations, two initial models IR and IC were prepared, 
respectively, by applying the following constraints: 
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FIG. 2: (color online) Comparison of S{Q) of Li20-Si02 RMC 
models based on crystalline (black dashed line) and random 
(blue dotted line) starting configurations with that of a Li2 0- 
Si02 MD system (red solid line). Neutron diffraction data is 
shown in panel (a) and X-ray diffraction data in panel (b). 



(i) Si is coordinated fourfold with O using a minimal 
neighbor distance of 1.4 A and a maximum neigh- 
bor distance of 1.8 A, This corresponds to a 100% 
fraction /si,o(4) = 1 of fourfold coordinated Si. 

(ii) The relative numbers of bridging (two Si neigh- 
bors) and non-bridging oxygens (one Si neighbor) 
is 37% and 60%, respectively. This corresponds to 
/o,si(2) = 0.37 and /o,si(l) = 0.60. 



(iii) Minimal atomic distances r„ 
are required. 



given in Table IIIII 



Intra-tetrahedral 0-Si-O angles and inter-tetrahedral Si- 
0-Si angles were allowed to evolve freely. The prepara- 
tion was run until the constraints (i)-(iii) were satisfied 
for at least 95% of the atoms (for given uncertainty pa- 
rameters, see below). 

After creating the initial models, the final refinement 
is done in order to obtain the best possible agreement be- 



TABLE III: Minimal atomic approach distances as used in 
the RMC modeling. 



TABLE IV: Weighting factors ac, used for the input quantities 
in the RMC modeling. 



Pair 


Si-Si 


Si-O 


Si-Li 


0-0 


0-Li 


Li-Li 


rmin in A 


2.8 


1.4 


2.5 


2.3 


1.7 


2.2 



quantity 
a 



S{Q) G^"^(r) /si,o(4) /o.si(2) /o,Si(l) 



10" 



10" 



10" 



10" 



10 



-3,5 



4 



tween the computed S{Q) and G^^^{r) from the RMC 
model and the calculated data from the MD-simulation. 
Both the real space as well as the reciprocal space data 
was used, since strong low-Q features in S{Q) empha- 
size the medium range order, while G^^^(r) shows well 
defined low-r features which emphasize the short range 
atomic order. The same constraints as in the prepara- 
tion of the initial models were applied also during the 
final RMC modeling. 

In the modeling, the input quantities, i.e. the total ra- 
dial distribution function G^^^{r), the total structure 
factor S{Q), and the fractions /si,o(4), /o,Si(2), and 
/o,Si(l) of differently coordinated silicon and oxygen ions 
are taken into account by an effective Hamiltonian of 
type Hcs = J2aXa/'^l> where Xa is the (unweighted) 
square deviation between the computed and the mea- 
sured (in the case of constraints presumed) value of the 
input quantity a. The weighting factors a a are sum- 
marized in Table IIVI For a detailed description of the 
algorithm we refer to the manual of the RMCA and and 
RMC-I--I- package which can be downloaded at [2^. 

In total five RMC models were produced. Two models 
were generated based on X-ray diffraction data using the 
random and the crystalline starting configuration (XR 
model and XC model) and another two are based on the 
neutron diffraction data (NR model and NC model). A 
fifth model was generated using both the X-ray and the 
neutron diffraction data starting from the random con- 
figuration (NXR). 

In Figs. [1] and the corresponding total structure fac- 
tors and total distribution functions are compared to the 
ones calculated from the MD structure. As expected, a 
good agreement is achieved through the RMC modeling. 



IV. COMPARISON OF STRUCTURAL 
PROPERTIES 
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A. Partial radial distribution functions 

In Fig. [3] partial radial distribution functions gij (r) (see 
Eq. [2]) are shown. Generally there is a reasonable agree- 
ment of the RMC models with the MD structure for Si-Si, 
Si-0, and Li-0. For goo the second peak is lacking or too 
weakly pronounced in the RMC models. This difference 
can lead to significant deviations of the RMC structures 
and the MD structure in the medium-range order regime 
(see also the discussion in Sec. lIVDp . While there is no 
dependence upon the starting configuration (crystalline 
or random), some differences can be seen between the 
RMC models based on the neutron diffraction data and 
those based on the X-ray diffraction data. These are 
most pronounced in gsiu and guu- The RMC models 
for which both the neutron and the diffraction data was 
used show the best agreement with the MD data. 

To summarize the result of the comparison of all RMC 



FIG. 3: (color online) Comparison of gij{r) of Li20-Si02 
RMC models based on crystalline (black dashed line) and 
random (blue dotted line) starting configurations with that 
of a Li2 0-Si02 MD system (red solid line). The RMC model 
based on neutron diffraction data is shown with a vertical off- 
set of 2, while the RMC model based on both neutron and 
X-ray difi'raction data is shown with an ofi'set of 4. In order 
to show the full first peak, the curves are scaled by the given 
factor on the left side of the vertical dashed line in each plot. 

models and the MD model, we calculated the integral 

rfS = / l5S -ffrid^ (6) 

over the difference between the partial radial distribu- 
tion functions of the RMC model a and the MD sys- 
tem. In the numerical calculation we integrated up to 
r = 10, which amounts to an integration to infinity, since 
c," ^ ^ 1 for r > 10. The resuhs shown in Table |V] 



TABLE V; Integrated differences of partial radial distri- 
bution functions, cf. Eq. (|6]). 



Pair 


AL; 


vr? 
AH 




INH 


1\TVD 

IN AH 


iU 


IH 




Si-Si 


LIO 


1.24 


1.42 


1.50 


1.21 


1.56 


1.56 


9.76 


Si-0 


L25 


1.29 


1.44 


1.44 


1.28 


3.58 


3.47 


6.58 


Si-Li 


1.47 


1.52 


2.36 


2.27 


1.53 


1.63 


1.83 


9.84 


O-O 


1.44 


1.47 


1.03 


0.99 


0.88 


2.29 


2.21 


5.93 


Li-Li 


2.03 


2.01 


1.05 


0.97 


1.14 


1.51 


1.58 


9.75 


Li-0 


1.46 


1.42 


1.24 


1.12 


0.91 


2.01 


2.04 


6.71 


Average 


1.46 


1.49 


1.42 


1.38 


1.16 


2.10 


2.12 


8.10 



allows us to quantify the quality of the RMC models 
relative to each other. It is surprising that on average 
the IC and IR models that do not involve information 
from scattering functions, are not much worse than the 
RMC models (XC, XR, NC, NR) with only one scatter- 
ing probe (X-ray or neutron). On the other hand, when 
information from both scattering probes is taken into ac- 
count, one obtains a significant improvement on average. 
We note that this improvement is associated with the 
fact that the NXR model exhibits good agreement for all 
individual partial radial distribution functions. By con- 
trast, when the analysis is based on one scattering probe 
only, quite large deviation can occur for certain partial 
radial distribution functions (see, for example, dgp^^^ and 
r/^C ^ 



B. Coordination numbers 
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FIG. 4; (color online) Comparison of histograms of coordina- 
tion numbers of Li20-Si02 for RMC models based on a crys- 
talline starting configuration using neutron diffraction data 
(black solid bars) and X-ray diffraction data (blue striped 
bars), and the original MD structure (red solid bars). The 
first atom type is the center atom, and the given distances 
are the radii of the coordination sphere. 



The first neighbor shell coordination number distribu- 
tions fij{N) are defined here as the fractions of atoms 
of type i which have N atoms of type j within a first 
neighbor shell radius r^ . Note that being neighbors in 
the sense of this analysis is not associated with having 
a chemical bond. To quantify the quality of the various 
RMC models a, the overlap 

OG 

was calculated. A value (jjf^ = 1 means perfect overlap 
between the distributions of the coordination numbers 
/g(A^) and n^°{N) of the RMC model a and the MD 
structure. 

The results are summarized in Table IVTl No significant 
differences in the 0^- are found between the RMC mod- 
els based on the crystalline starting configuration and 
the random starting configuration (with one exception 
for Si-Li, where a larger difference is observed between 
the RMC-XC and RMC-XR models). This suggests that 
the quality of reproducing coordination numbers is inde- 
pendent of the starting configuration. A significant dif- 
ference between the neutron based and the X-ray based 
RMC models is found in the Li-0 coordination numbers. 



where the better overlap for the neutron based model can 
be traced back to the higher relative weight of lithium in 
the neutron diffraction functions. An improved overall 
agreement is achieved when using both X-ray and neu- 
tron diffraction data, though the two most significant 
discrepancies (Si-Si and 0-0) are still there. It is also 
informative to take a look at the overlap numbers of the 
initial models that are on the constraints only. These 
are comparable in quality with the RMC models, which 
in addition take into account the information from one 
scattering probe. As for the partial radial distribution 
functions discussed in the previous Sec. IIV Al the RMC- 
NXR shows a clear improvement compared to the initial 
RMC models IC and IR. 

Figure 2] shows a detailed comparison of the coordina- 
tion number distributions of the RMC-NC and RMC-XC 
models with the MD model. For the distribution /sio(^) 
not shown in Fig. [4l we obtained a very good agreement 
which essentially results from the constraint that silicon 
atoms must have 4-fold coordination. The most strik- 
ing discrepancies between the RMC models and the MD 
structure are found in foo{N) and /siSi(-/V). The MD 
model shows a clear bimodal distribution with maxima 
at 3 and 6 neighbors (corresponding to non-bridging oxy- 
gens and bridging oxygens) in tiqq{N), while the RMC 
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TABLE VI: Overlap of coordination numbers in percent (see 
Eq. [7]); first neighbor shell radii are given in the second col- 
umn. 





[A] XC XR NC NR NXR IC 


IR 


C 




9 n 




yo.o 


QQ 7 


QS A 
yo.^ 




99.7 98.8 


93.2 




o n 
2.U 


yy.o 


99.6 


99.0 


99.6 


99.1 


99.6 


99.6 


91.2 




2.8 


69.8 


68.2 


83.0 


85.1 


87.7 


69.3 


69.8 


48.8 


O-Li 


2.8 


84.5 


82.5 


84.9 


84.9 


86.6 


83.1 


84.5 


57.6 


Si-Li 


3.8 


90.4 


80.9 


77.8 


78.5 


83.8 


80.6 


74.5 


47.7 


Li-Si 


3.8 


89.5 


87.9 


87.1 


86.7 


89.2 


84.7 


81.7 49.0 


Si-Si 


3.5 


64.0 


65.1 


64.5 


61.9 


62.3 


66.0 


64.1 


63.9 


O-O 


2.9 


52.9 


52.0 


54.2 


54.3 


58.9 


51.6 


51.1 


89.9 


Li-Li 


3.5 


82.4 


81.5 


79.5 


81.2 


92.8 


84.4 


84.9 


52.6 


Average 




8L3 


79.6 


81.1 


81.2 


84.0 


79.9 


78.8 


66.0 



models have a broad smooth distribution. On the other 
hand, ng^^{N) is much narrower than ng?gi(A^). These 
findings suggest that the short-range order of the RMC 
models corresponds quiet well to that of the MD struc- 
ture, but that the medium-range order, and particularly 
the structure of the Si-0 network, has significant differ- 
ences. We note, that there are virtually no differences 
between /siSi(Af) and foo{N) among the five RMC mod- 
els. 

In summary we can conclude that most features in the 
coordination number distribution are already captured 
by the constraints. This may not be surprising, since 
coordination numbers for Si and O have been used as in- 
put requirements together with the rather high density of 
the system. As a consequence, there is not much freedom 
for the coordination numbers between other types of ion 
pairs. 



C. Bond-angle distribution 

We calculated bond-angle distributions for intra- 
tetrahedral angles (0-Si-O) and inter-tetrahedral angles 
(Si-O-Si) and found that all RMC models have essentially 
the same bond-angle distributions. Differences lie within 
the statistical spread. 

In Fig.[5]the distributions for the RMC-NR and RMC- 
XC model are compared to that of the MD structure. 
The Si-O-Si bond-angle distribution of the RMC mod- 
els agrees well with the MD data. The intra-tetrahedral 
bond-angles, on the other hand, are much broader dis- 
tributed in the RMC models than in the MD struc- 
ture. This impression can be quantified by calculating 
the mean angles a(Si-O-Si) and Q:(0-Si-0) as well as the 
standard deviations Aa(Si-O-Si) and Aa(O-Si-O). It is 
found that the mean angles of all RMC models agree very 
well with the MD values, while the standard deviations 
are larger by a factor of two, see Table IVlTl This finding 
corresponds to the deviations observed in the partial ra- 
dial distribution functions in Fig. [3] and the coordination 



TABLE VII: Mean values and standard deviations of bond 
angles. 

MD XC XR NC NR NXR IC IR C 

a(Si-O-Si) 141.2 139.8 140.7 137.4 138.6 140.5 138.0 138.5 125.7 

Aa(Si-O-Si) 13.7 13.1 13.6 14.5 14.5 13.1 14.9 15.5 0.5 

a(O-Si-O) 108.8 108.7 108.7 108.7 108.8 109.0 108.2 108.2 108.9 

Aa(O-Si-O) 5.9 13.1 12.9 12.7 12.4 11.3 15.6 15.9 3.2 



number distribution of 0-0 in Fig. |4l There, distinctive 
features of the MD data, as the second peak in and 
the bimodal distribution in , are not well reproduced 
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FIG. 5: (color online) Bond-angle distributions of the MD 
structure (red solid lines) and RMC models based on a crys- 
talline starting configuration using neutron diffraction data 
(black dashed lines) as well as based on a random starting 
configuration using X-ray diffraction data (blue dotted lines). 
In the left panel the distributions of intra-tetrahedral angles 
(O-Si-O), and in the right panel the distributions of inter- 
tetrahedral angles (Si-O-Si) are shown. 



D. Ring-size distribution 

In order to compare the topography of the glass- 
network we determined ring-size distributions for each 
model. Here rings and their size are defined in the fol- 
lowing manner: 

(i) A Si-atom and an 0-atom are considered neighbors 
if their distance is smaller than 2.0 A (using closest 
image convention). 

(ii) For each Si-atom 7, the smallest closed loop of al- 
ternating neighboring silicon and oxygen atoms is 
determined, which entails the Si-atom 7. 

(iii) The size of the ring is equal to the number of its 
Si-atoms. 
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FIG. 6; (color online) Comparison of ring-size distribution 
of Li20-Si02 RMC models based on a crystalline starting 
configuration using neutron diffraction data (black wide bars) 
and based on a random starting configuration (blue striped 
bars) with that of a Li20-Si02 MD system (red narrow bars). 



FIG. 7: (color online) Picture of the RMC model based on a 
crystalline starting configuration. Lithium atoms are marked 
as small yellow spheres, Silicon atoms as medium sized green 
spheres, and Oxygen atoms as large red spheres. In the center 
a linear chain of neighboring Si and O is marked in dark blue, 
which is a leftover of the starting configuration. 



The maximum number of rings equals the total number 
of Si-atoms. However the number of rings is generally 
smaller, since there are a number of Si-atoms for which 
no ring is found (e. g., for an isolated Si04 tetrahedra), 
and two different Si-atoms can be associated with the 
same ring. 

In Fig. [HI the ring-size distributions of the MD system, 
the RMC-NC and the RMC-NR-model are shown. No 
data is shown for the X-ray and combined data based 
RMC models, since their ring-size distributions are prac- 
tically the same as for the neutron diffraction data based 
RMC models. Indeed the ring-size distributions almost 
do not change compared to those of the initial models. 
On the other hand, there is a clear dependence upon the 
starting configurations. While most rings (40%) of the 
RMC-NR model are of size three and four, the RMC- 
NC model has a high number of rings of size 4 and 6. 
The latter is more in line with what is found in the MD 
system. 

The RMC-NC model has a high number of rings of a 
size larger than 10. Examining these large rings in more 
detail reveals that most of them are actually straight lin- 
ear chains penetrating the system parallel to one system 
axis (they are seen as rings due to the periodic bound- 
ary conditions, and can be found in both the RMC-NC 
and the RMC-XC models; see also Fig. [T]). Such straight 
chains are not seen in the RMC-NR and the RMC-XR 
models and in the MD structure. There are also large 
rings and chains in these models, but those are generally 
much more twisted than in the RMC models based on 
the crystalline starting configurations. 



V. SUMMARY AND CONCLUSIONS 

The RMC method successfully reproduces many 
salient features of the local structure of the original MD 
system. Some differences are found in the partial ra- 
dial distribution function of 0-0, and in the coordina- 
tion number distributions of Si-Si and 0-0. With respect 
to the structure beyond nearest neighbor distances the 
RMC models are less predictive and therefore cannot be 
expected to capture the medium range order properly. 

Comparing RMC models based on X-ray and neutron 
scattering data revealed no significant differences. More- 
over, we found that the additional consideration of scat- 
tering data from one type of probe (either X-ray or neu- 
tron) gives only a modest improvement over the initial 
RMC models that are based on geometric constraints 
only (number density, minimal pair distances, some coor- 
dination numbers). The situation becomes much better, 
however, when including the information from both types 
of scattering probes. 

Most structural properties of the RMC models do not 
depend sensitively on the starting configuration (crys- 
talline or random). Even the ring-size distributions do 
not differ that much. However, taking a closer look at 
the rings, revealed that the RMC models based on the 
crystalline starting configuration exhibit straight linear 
chains penetrating the system. These straight chains 
are remnants of the crystalline starting configuration and 
their occurrence is not reflected in the other structural 
properties studied. In particular, there are no differences 
in GP°P(r) and S{Q) between the RMC models based 
on the crystalline and the random starting configuration. 
These findings show that one should check carefully if a 
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feature of interest in RMC models is only a product of (DFG Grant number MA 1636/3-1) and by the NSF 
a particular starting configuration or if it can be repro- (NSF DMR Grant number 0710564). 
duced using totally different starting configurations. 
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